Wilson-KadanofF Renormalization Group in Higher Orders: 
One-Dimensional g-ology Model as an Example 

Gennady Y. ChitovQ 

Centre de recherche sur les proprietes electoniques de materiaux avances, 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada JIK 2R1 and 

Physics Department, McCill University, 
Montreal, Quebec, Canada H3A 2T8 

Claude Bourbonnais 

Centre de recherche sur les proprietes electoniques de materiaux avances, 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada JIK 2R1 
(Dated: February 1, 2008) 

Abstract 

We apply the standard Wilson-Kadanoff (WK) momentum-space Renormalization Group (RG) 
scheme for the g-ology model of one-dimensional fermions. By explicitly carrying out calculations 
at the two-loop level, we show how the RG flow equations can be derived from the summation of the 
cascades of contractions generated by the effective action's mode elimination at each infinitesimal 
step of the WK procedure. The rules for selecting these series of cascades appear naturally as 
a consequence of the WK scheme "on-shell" kinematic constraints and conservation laws. The 
relation between the present RG approach and the field-theoretic schemes used in earlier related 
studies is analysed. Generalizations for other models and/or higher dimensions are formulated. 

PACS numbers: 05.30.Fk, Tl.lOPm, ll.lO.Hi, 05.10.Cc 



*Present address: Department 7.1-Theoretical Physics, University of the Saarland, Saarbriicken 66041, 
Germany 



I. INTRODUCTION 



The Renormalization Group (RG) approach, in its most enhghtening formulation due to 
Wilson m, 0], is the theory designed to handle fields (quantum or classical) fiuctuating over 
range of momentum (energy) scales.^ In the momentum-space formulation, this range of 
scales where fiuctuations matter is bounded by some initial ultraviolet (UV) cutoff Aq (not 
necessarily finite) provided by the problem under consideration, and Aq is the scale where 
the RG procedure of the successive mode integration starts. For the most problems where 
RG was successfully applied (e.g., phase transition, field theories) the low-energy sector of 
interest resides near one point of the momentum space which can always be mapped to the 
origin of that space. 

In practical loopwise calculations, the original WK infinitesimal scheme seems however 
to appear too involved beyond one- loop level, and, e.g., for the Wilson's original problem of 
phase transition in higher orders various field-theoretic RG approaches were applied [Q. In 
field theory the cutoff Aq is essentially considered as a "troublesome" parameter to get rid 
of, either by eventually taking the limit Aq — > oo, or, e.g, by using another regularization 
scheme, like the dimensional one where the cutoff is set to infinity from the start. 

Contrary to the WK scheme, which specifically stipulates that at each RG step the 
momenta to be integrated out lie within infinitesimally narrow shell in the momentum 
space, in the field-theoretic RG only overall momentum conservation is taken care of, while 
separate momenta may lie outside of the (hyper)sphere of radius A (centered at the origin). 
This situation is well known, and contributions from the states lying outside of the action's 
phase space at the distance of several radii A- or, in other terms, violating phase space 
constraints by (9(A)- do not affect the leading terms of the RG flow. 

However, in condensed matter problems the Wilsonian effective action's scale Aq must 
be taken seriously. A priori it cannot be set to inflnity or even to the largest scale provided 
by the underlying microscopies, i.e., inverse lattice, interatomic, etc, spacing. Even for 
a problem with the low-energy sector localized near the origin of the momentum space, 
aforementioned violations of the phase space constraints result in double counting of the 
degrees of freedom lying outside of the Wilsonian phase space (i.e., the degrees of freedom 
which are already "integrated out" in some way before we reached the Wilsonian scale 
from the microscopic scale). This can affect the next-to-leading terms, as well as flnite 
renormalizations of the action's parameters. 

For the important problem of fermions having a Fermi surface - that is a low-energy 
sector residing near a manifold in the momentum space - accounting for the modes within 
the phase space only during the RG procedure becomes a rather non-trivial task, as was 
explained by Shankar 0. Moreover, the error of violating the phase space constraints can 
be of order 0{Kf), where Kp is some characteristic size of the Fermi surface, even if the 
low-energy modes we are integrating out reside within a narrow shell A <^ Kp. 

Note that the answers for many of the most interesting problems of the condensed-matter 
fermions lie in calculating the abovementioned "next-to-leading terms" , where the differences 
between various RG approaches matter. For instance, the problem of the Fermi-Liquid vs 
non-Fermi-Liquid regimes of the normal state of the cuprates, from an RG point of view is 



Since the Wilson RG theory was strongly motivated by Kadanoff's more intuitive original approach of 
successive averaging of the spin Hamiltonian in real space ||^ , it is often called the Wilson-Kadanoff (WK) 
RG, the term we apply in this study. 



2 



the difference in the behavior of the irrelevant (self-energy) terms in both regimes. So the 
choice of the optimal RG scheme and the control over the subleading terms are the issues 
far from being academical. 

The main claim we want to make from the considerations of this paper is that for a 
condensed-matter problem cast in terms of the Wilsonian effective action with a physically 
meaningful UV cutoff Aq, the WK RG provides the only scheme definitely reliable in all 
orders of the loopwise expansions, where all modes from the phase space are integrated 
once and only once. (Strictly speaking, there are no new results in the RG equations we 
derive below, since for the ^f-ology model they have been known for a long time, either from 
field-theoretical RG scheme H, |^ or more recent WK RG formulations p). We address 
the issues of how to practically employ the WK RG in higher orders, to which technically 
simpler scheme the WK RG can be reduced and to what extend, by studying the g-ologj 
model of ID fermions. 

The model we consider is the simplest case of the condensed-matter fermions where the 
Fermi surface manifold is just a set of two points. By going deeper into a recent suggestion 
to this problem 0, we apply a WK RG transformation assuring that at each step of the 
procedure all momenta lie within the infinitesimal shell of integration and satisfy the con- 
servation laws, so there are no contributions from the extra states nor double counting. We 
show how the higher order renormalizations of the effective action's parameters appear from 
the summation of the cascades of contractions generated at each step of the WK procedure. 
We don't go beyond the two-loop level, which is apparently the optimal one for feasible cal- 
culations. At this level, we show that the summations over the cascades in the WK scheme 
is equivalent (in the fixed-point limit of the decreasing cutoff A — 0) to applying the cutoff 
regularization for the field-theoretic RG with the extra states excluded from the phase space. 
This result, to the best of our knowledge, was not explicitly shown before. Considerable 
simplifications for our analysis came from introducing the chiral matrix formalism. This is a 
naive (for the simple ID case) first attempt to build into the technique the disentanglement 
of geometry (characterized by the "large" rigid scale Kp) and the low-energy sector per se 
(characterized by the decreasing scale A), and to map the latter to the same origin. This 
mapping makes the exclusion of the extra states particularly simple. 

The rest of this paper is organized as follows. Sections [II A\ - fl D] are introductory. Section 



II A| defines the effective action. In Section III B| the chiral matrix formalism is introduced. In 



Section [II D| we explain the terms and notations of the Wilson-Kadanoff RG as it is applied 
in this paper. In Sections [III A[ - pTB[ we explain the procedure and give the results for the 
two-loop self-energy calculations, and in Sections [III C[ we establish the relationship between 
the field-theoretic RG and the present approach. Section IV contains the two-loop results for 
the couplings flow equations. The conclusions and generalizations for other models and/or 
higher dimensions are presented in the final Section 0. 



II. MODEL AND FORMALISM 



A. SU (A^)- invariant d-dimensional ^/'^-action 



We treat d- dimensional interacting fermions at finite temperature in the standard path 
integral Grassmannian formalism |T^. The 5't/(A^)-invariant fermionic ^/'^-action along with 



the related quantities and notations of this subsection have been introduced earlier [|TT], and 
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we recall these here in order to make the paper self-contained. The partition function is 

Z = j VtpV'ip e^o+^'"' (1) 

The free part of the action is 

^0=/ V^,(l)K + /i-e(Ki)]^,(l) (2) 
7(1) 

and {d + l)-dimensional vectors and integrals mean 

... " ^/(^? 

(i) = {K,,uj,) (3b) 

where f5 is the inverse temperature, /i the chemical potential, uji the fermion Matsubara 
frequencies and ipa{i) a A^-component Grassmann field with a spin (flavor, if 7^ 2) index 
a. Summation over repeated indices is implicit throughout this paper. We set fc^ = 1 and 
h = 1. The S't/(A^)-invariant quartic interaction is 

Sint = -] f ^„(l)^;3(2)^,(3)^,(4)C/^f (l,2;3,4)5('^+i)(l + 2-3-4) (4) 

4 7(1,2,3,4) 

Here the conservation of energy and momentum is enforced by the symbolic delta function 
+ 2 - 3 - 4) = /3(27r)'^5(Ki + K2 - K3 - K^)6^,+^,^^,^^,,o (5) 

Note that for the lattice fermions the delta-function above conserves momenta up to an 
inverse lattice vector. The physical case of the spin-| electrons is recovered by setting 
N = 2. We decompose the interaction by factorizing its symmetric and antisymmetric 
frequency-momentum- and spin- dependent parts as 

= f/^/^f + U'T^f (6) 

where the functions U^^^ are symmetric/antisymmetric under exchange of the left or right 
pairs of their variables. Two operators / and T, which are respectively symmetric and 
antisymmetric in the spin space, are defined as follows: 

/"f = daeSfSy + Sa-ySfSe (7a) 
T^/ = ^aeh-y - ^'^ihe (Tb) 

The one-particle-irreducible (IPI) four-point vertex F, defined in the standard way, com- 
prises two components (F"^, F*^) and can be written in the same manner as (^. Along 
with such representation showing explicitly its antisymmetry, there is another one which 
separates the vertex into charge and spin components via 

f ^ F^f = F^J-f + F^T^f = -^T^5^,8,, - \t^\1^\% (8) 
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The latter representation is often more convenient since charge and spin components decou- 
ple in many practical calculations. A"^ (a = 1, A^^ — 1) in Eq.(^ are Hermitian traceless 
generators of the SU{N) group, coinciding with the Pauli matrices for N = 2, and normal- 
ized such that 

tr(r A^) = 25"* (9) 
They also satisfy the following Fierz identity: 

The components of the vertex in different representations are related as 

= (A^ - l)r^ - (A^ + l)r^ (11a) 

r, = -r^-r^ (iib) 



B. Chiral matrix formalism for the ID effective action 

The construction of the low-energy effective fermionic action as a simplified form of the 
generic ■^/''^-action (0,^) reduces in ID essentially to the three key steps: extraction of a 
finite set of (marginal) couplings; decomposition of the physical electron ip-field into two 
chiral fields "near" left and right Fermi points; linearization of the one-particle spectrum 
e{K) around these points. There exists an extensive literature on the derivation of the ID 
effective action models (Tomonaga-Luttinger, g-ologj), and on many subtleties it involves 
(see, e.g., Ref. and more references therein). So below we simply present notations, 



approximations and highlight the main points. The chiral matrix representation we derive 
here is in fact only a way to re-write the well-known gf-ology model. (For details see the 
next subsection.) However this matrix formalism and related to it diagrammatics help to 
facilitate bookkeeping of diagrams and handling their contributions to different couplings, 
and, eventually, to generalize the g-o\ogy results for other cases. 
We decompose ID momenta as 

= p,Kf + ki (12) 

where Pi = sign{Ki) = ±1 corresponds to the right (R) /left (L) Fermi points. In the 
quartic interaction function U{1, 2; 3, 4) (^) we retain only its dependence on the momenta 
at the Fermi points, neglecting that on fcj-s and frequencies . The SU (N) and exchange 
symmetry constraints leave us with four independent marginal couplings: 

= U'^{LR;LR) = -U'^{RL;LR) = -U^{LR;RL) (13a) 

/ = U^{LR;LR)= U^{RL;LR)= U^{LR;RL) (13b) 

c/4 = U^LL; LL) = U^RR; RR) (13c) 

g3 = U^{LL;RR) = U^{RR;LL) (13d) 

Note that the Umklapp coupling g^ is operative only in the half-filled case when AKp = 2tx 
(we set the lattice constant a = 1). For the formalism we develop below this case is considered 
as the most general one, while the non-commensurate filling can be recovered by setting 
^73 = 0. 
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We introduce left (right) chiral Grassmannian fields as: 

MtKf + ki, uji) = ^^'^{ki, Ui) (14) 

Small vectors 

fciG[-Ao,Ao] (15) 

are restricted by the bare momentum cutoff Aq of the effective action. Then the effective 
quartic interaction [cf. Eq.(^] can be represented as 

S':^ = -\l 4^5(^+^^(1 + 2-3-4) (16) 

^ ^(1,2,3,4) 

where the transition from bold to thin variables corresponds to h-* k^ [cf. notations 
(IgpD], and 

4J = ^^(l)#(2)^,^(3)^,^(4)[ ■ J^f + g, ■ T^f ] + {L ^ R) 
+ V^^(l)V^^(2)V^f (3)^f (4)[ ■ J^f + g, ■ T^f ] + {L ^ R) 
+ ^^(l)^J(2)^^(3)^f (4)[ g^ ■ J^f + g' ■ T^f ] + {L ^ R) 
+ ^^{l)^^i2)i;^m^iA)[-g^ ■ J^f + / ■ T°/] + (L ^ R) (17) 

Note that the delta- function in (|1^), contrary to its "bold" counterpart (^, conserves the 
small momenta ki exactly, since the lattice effects at half-filling (^fa-terms) are explicitly 
taken into account in (p!?]). 

Let us now introduce four 2x2 matrices f such that is the unit matrix, and for 
i = 1,2,3 correspond to the Pauli matrices. Such distinct notation helps to avoid confusion, 
since f operates in the chiral, and not in the spin space. We attribute an extra chiral index 
to the Grassmannian field as L i— 1, i? 2. Then four terms in the last two lines of Eq.(|T^ 
can be written as: 

= ■ /,f ■ tl,,t'^'s^ + / ■ T^f ■ ti,,,t'y,,] C'(l)^f (2X(3)Vr (4) (18) 
while the first two lines in that equation give 

= T^f ■ [g, ■ t%>t'y,, + g^m ■ tl,p,t'^,,,] C'(l)V^?'(2)<(3)V^|'(4) (19) 

where ^ 

9p/m = -{g4 ± gs) (20) 

Thus 

4ff — -^fb + -^uc (21) 

The one-particle Green's function is spin-independent in the model we consider, but it 
depends on the chiral index: 

-(V^^(l)V^f (2)) = G«(l)<5.^5,,<5(i+i)(l - 2) (22) 

The linearized bare Green's function is given by 

Go^^ikn, UJn) =iuJn + fl- (^{pc^Kp + /Cn) ~ i^n - PaVpkn (23) 
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where 

= 2(a - 1) - 1 (24) 

Now we can define the IPI four-point vertex f in this formahsm. We will be interested in a 
particular limiting form of this vertex when its dependence on the external momenta lying 
at the Fermi points (parametrized by the chiral indices) is retained only, while (fc„, a;„) are 
discarded. In the lowest order f can be read off from £gg (|T^,|r^,^). As follows from the 
group-theory arguments, this structure of the vertex is preserved by interaction C^^: the 
latter results in only the renormalization of couplings, but no new couplings are generated. 
Thus, the following operator expansion with four couplings holds non-perturbatively: 

f = -g^ / ® P ® F + f ® [g^ f f + c/p ® t° + gm P ® (25) 



In order to keep the notations compact, we do not change them in (|25|), however couplings 
are renormalized by interaction comparatively to their bare values (|T3|). To the same end 
we will not indicate explicitly the chiral and spin indices from now on (unless there is a 
possibility of confusion), incorporating them into the "thin variables". So we imply 

Cn. + C,.c ^ [cf. Eqs.(|l|,[l9|)] ^ r(o)(12; 34)V^(1)^(2)V^(3)V^(4) (26) 

Comparison of Eqs.([T3JT^123) with the operator form (P5|) clarifies how we attribute the 



chiral and spin indices for the matrix elements of T. Note the explicit antisymmetry of the 
vertex r(12; 34) under exchange of the left or right pairs of variables in the representation 

(El)- 

We can recover couplings from the vertex by appropriate convolutions, i.e., 

-^(4-r(12;34)-4^-4) (27a) 
l(t°,-r(12;34).T,^,3-4) (27b) 



9p = 

and other two couplings are obtained from Eq.([27|b) by obvious substitutions, namely. 



C. Contact with the (^-ology model 

Here we want to make a contact between the chiral matrix formalism and more commonly 
used g-o\ogy language (see, e.g., Ref. [^]). No surprise that this subsection's formal manip- 
ulations with the fully (anti)symmetrized "chiral" action amount essentially to reducing it 
to the non-symmetrized form. The latter is usually used as a departure point of various 
g-o\ogy models. 

Let us show first how the charge- and spin-operator interaction action, involving forward 
and backscattering couplings called g2, gi respectively in the g-ologj nomenclature, can be 
recovered from the part £fb (|TB|) of the effective quartic interaction (PT|). £fb can be split in 
two terms by the following projection in the chiral space: 

^fb = {[l-t%t%]+t%t%){-g^itl,f^, + g'ftl^^^^^^ 

= -^csl + ^cs2 (28) 



7 



Above we indicate explicitly the chiral indices only, and [V'^Itv stands for the product of four 
■j/'-fields "normally-ordered" as on the r.h.s. of Eg. (p!8D . Using the Fierz identity (10) for the 
SU{2) chiral operators f, one can prove that 



-StlAs + tl,Ae)i-9^I + g'fmN (29) 



1 

= 2' 

It can be also shown that the above expression reduces to 

'(^csi = tipi-g^i + (7^T)r (l)V'^(2)V'''(3)r (4) (30) 

In the same manner we find 

Ccs2 = tipig^i + /f )r (l)V^^(2)r (3)V^^(4) (31) 

Using the symmetry properties of the spin operators and changing the dummy variables, it 
is easy to see that the term £csi results in the same expression for the action [cf. Eq.(|l^)] 
as £cs2, so Sfh = 2S'cs2- Applying then the Fierz identity (|1^) [cf. also definitions (0)] for the 
operators I,T in case of SU{2), and using the conservation law for momentum-frequency 
such that 3 = 1 + Q, 4 = 2 — Q, we get 

-^fb = —tl^ I {gcP''{Q)p^{-Q) + Ag,S:{Q)SP,{-Q)) (32) 

J Q 

where we introduced the operators of the chiral (i.e., left/right) charge and of the chiral 
a-component of spin as 



p^iQ) = Jr^ilWl + Q) (33a) 

Note that in the above formulas there are no sums over a. We have also defined the charge 
and spin couplings: 

gc= g'- Sg^ (34a) 

gs = -g' - g^ (34b) 

It is worth noting that although there is no one-to-one correspondence between ID cou- 
plings (|TB|a,b) and the vertex components F"^, F'^ [the latter in ID are rather products of 
corresponding couplings and operators in the chiral space, as can be seen from comparison 
of Eqs.(|8|, |25|) ], the definitions for the components which couple separated charge or spin 
modes ([TlD hold for the chiral ones (|3^ . It is also useful to give the relationship between 
our constants and the standard (dimensionful) (^-ology couplings: Qc = 2g2 — gi, gs = —gi- 
The term S^c of the total quatric interaction generated by the contact coupling at the 
same Fermi point (5^4) and the Umklapp coupling ((73) can be written as 



^uc= - ^^(p"(Q)p"(-g)-45:(g)5:(-g)) 

fap I {^.J^e - \%K)r,ii)€iX + Q) ■ C(2)€(2 - Q) (35) 



^3 

4-2' 
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The way we write the Umklapp contribution is to reinstate the known fact that since 
couples fields with different chiralities, it cannot be represented in terms of the local chiral 
charge or spin operators. 



D. Wilson-Kadanoff renormalization group for the chiral effective action 

In this study we apply the original Wilson-Kadanoff momentum-space Renormalization 
Group scheme [|I|, |^ to the effective action with the quartic interaction containing four 



coupling constants. This interaction in the chiral formalism is given by Eqs. ( p5| , p6D . In 
that formalism all the constraints coming from the original two-point geometry of the Fermi 
surface, and consequently, occurrence of the two types (left/right) of fermion fields and of 
Green's functions, are taken care of by the appropriate summations over chiral indices. In 
applying the RG procedure we successively integrate out the Grassmannian fields in the 
phase space of the small momenta ki (|12|). The latter, like in the quantum or classical (p'^- 
theories, can be seen now as having common origin towards which we eliminate modes in 
the /c-space. 

We start with the effective action ( P3| , pBD having momenta restricted by the initial (bare) 
ultra-violet (UV) cutoff Aq (|TH). After steps of the mode elimination {N ^ 1), we end 



up with the cutoff A chosen such that: 

^»l (36) 

(We choose the RG scheme without rescaling of the cutoff.) Each step consists of mode 
integration within the infinitesimal shell A in the A;-space, where 

A ^ (37) 

(the equidistant step is chosen for simphcity only) such that after n steps the UV cutoff A„ 
is 

A„ = Ao - riA, \<n<N {An = A) (38) 
We define the n-th shell in the momentum space Sfi as 

Sn = [-A„_i, -A„] U [A„, A„_i] = s;; U 5+ (39) 

In the final results given in the following sections we will smoothly take continuous limit 
A — >■ {N — >■ oo) such that relationship (^) holds. Then the lowering cutoff A will be 
determined by the continuous RG flow parameter / deflned as: 

A(/) = Aoe-' (40) 

As it was discussed in the closely related earlier study of the g-o\ogy model 0, the straight- 
forward application of the Wilson-Kadanoff RG beyond one-loop level encounters some diffi- 
culties. They stem from the constraint on the momenta to be integrated out at the n-th step 
to lie within the shell iS„, and on the total momentum fiowing through a given graph to con- 
serve. It was shown in Ref. that the solution of this problem consistent with the logic of 
the WK RG can be found when multi-loop graphs (the standard ones for the field-theoretical 
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RG) renormalizing action's vertices are obtained through cascades of contractions at differ- 
ent steps of the RG procedure. In this study we refine the "cascade scheme" of the RG 
calculations put forward in Ref. for the ^f-ology model. We show that the aforementioned 
constraints on the momenta can be cast in terms of the "selection rules" for the numbers 
of steps entering a given cascade of contractions. Summations over all allowed types of 
cascades provides the results particularly illuminating for understanding the origin of the 
(logarithmic) equivalence of the WK and field-theoretic RG versions. 



III. SELF-ENERGY AT THE TWO-LOOP LEVEL 

A. Cascades of contractions and kinematic constraints 

We use the standard definition for the field renormalization constant relating it to the 
self-energy (assuming that the analytical continuation from Matsubara frequencies is done): 



,-1 ^ ^ _ dJ:aik,uj) 



It is easy to show that in our model Z does not dependent on the chirality index. The first 
correction to Z comes from the two-loop "sunrise" self-energy diagram. However, in the 
WK RG scheme such a diagram if obtained in the second order by a contraction of three 



legs of the four-leg interaction (|T^) with itself, gives a vanishingly small contribution. This 
is due to the momentum constraints discussed above. (See also Ref. ^ for more details.) 
The graphs we are looking for can come only from two types of cascades of contractions 
generated by the mode elimination of the original action (|T^) : 3-step and 2-step cascades 
depicted in Fig. 1. 3-step cascade (Fig. la) is generated (in the second order over interaction) 
as follows: (i) at the n-th step of the mode elimination the four-point vertex is contracted 
with itself to generate the six-leg vertex; (ii) m steps after [i.e. at the (n + m)-th step] 
two legs of this vertex are contracted; (iii) at the {n + m + k)-th step two other legs are 
contracted resulting in the final correction to the two- leg vertex (i.e., one-particle Green's 
function). 2-step cascade (Fig. lb) is generated as follows: (i) at the n-th step two pairs 
of legs of the four-point vertices are contracted; (ii) at the (n + m)-th step two more legs 
are contracted. The above description of the cascades was very schematic and illustrative. 
We did not take care of the kinematics of diagrams (that is why we did not indicate the 
momentum directions in the graphs of Fig. 1), and the integers numerating the steps are 
unrelated and free to take any values restricted only by the total number of steps of the RG 
procedure, i.e., l<n<n-\-m<n-\-m-\-k< N ? 

To proceed with the further analysis of the pertinent cascades let us consider the renor- 
malization of the one-particle Green's function at = A. As the RG flows towards the fixed 
point (A — s> 0) this quantity will provide us with the field renormalization constant Z (p|). 
For this purpose we consider the terms generated by the two types of cascades described 
above renormalizing the coefficient of the action's Grassmannian binomial il){l)ip{l). Dia- 
grammatically they are shown in Fig. 2. Analysis reveals that there are three kinematically 
non-equivalent graphs (a, 6, c in that figure) generated by the 3-step cascades and two graphs 



^ Note that 2-step cascade or "1-step cascade" (i.e., simple contraction) can be considered as special cases 
of a 3-step cascade if we allow the integers m, k to take zero values. 
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{d, e) from 2-step cascades. We choose the external momentum 

h = AN = A (42) 

and the other external variables (chirality, spin, frequency) are not specified at this point. 
For the further applications we approximate 

A = LA (43) 

where L is an integer (1 -C L <^ A^). The error 0{A) of the representation (^) disappears 
when we eventually take the limit A — 0. 

Let us start with the kinematics of the first three graphs. We find that for each of 
these graphs there are two solutions compatible with the constraints imposed by the WK 
procedure and the momentum conservation law. They are given in Table I. 

Table I: Solutions for the momenta of the three graphs in Fig. 2 compatible with the kinematic 
constraints. 



Graph 


a 


a 


b 


b 


C 


c 


Solution 


I 


II 


I 


II 


I 


II 


m'= 


m 


m-2L 


m 


m-2L 


m 


m-2L 


k2 e 




'-^N-m+2L 






c- 


c+ 


ks e 








^N-m+2L 




^N-m+2L 


k^ e 


C+ 


c- 


C+ 


c- 







The ranges for the integers n, m satisfying the solutions I, II are: 



I {a,b,c) : 1 <n<N -2m; I < m < — (44a) 
II{a,b,c): 1 <n< N + 2{L-m); 2L <m< — + L (44b) 

Thus in order to obtain the total contribution of these cascade-generated graphs to the 
self-energy we have to integrate two independent momenta in each graph over available range 
(i.e. two possible configurations I, II for each graph) and to sum over all allowed types of 
cascades parameterized by {n,m) within the ranges (^If). (And of course, there are always 
straightforward summations over other variables, i.e., chirality, spin, frequency which take 
all available values.) 

We have not shown in Table 1 available kinematic solutions for the 2-step cascade graphs 
d, e, since their contribution is negligible in comparison to that of a-c. Indeed, as can be 
shown from the results below, the graphs d, e as well as a-c at a fixed m give the contributions 
0{A/A) to the constant Z. However the extra sum over m renders the latter logarithmic 
O(ln^), while the former do not have this extra sum and thus strictly disappear in the 
limit A — s> 0. For the same reason we drop some extra solutions for the graphs a-c which 
occur only under special relationships between n and m (i.e., there is no double sum over 
n,m): they give vanishing contributions.^ 



^ In the light of the previous comment (see the footnote on p. ^ it is clear that special configurations (or 
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B. Two-loop calculations 



Now we can proceed with the RG calculations via a straightforward evaluation of the 
graphs from Fig. 2. For certainty we fix the external chiral index of those graphs to be 
2, i.e. corresponding to the right component of the Grassmannian field. Since graph's 
contributions in the spin space are diagonal and spin-independent, the external spin can be 
"up" or "down". After summation over chiral and spin indices we obtain for the sunrise 
graph contribution F_|_(l) (here plus stands to indicate our choice of the external chiral 
index) : 

F+(l) = / \\{gl + 3^7,') ■ + 2gl . ri + 2gl ■ ^^1 b^^'^^\\ + 2-3-4) (45) 

-'(2,3,4) '-^ ^ 

where thin variables signify momenta and frequencies, i.e., the only ones left, and 

= G_(2)[G'_(3)G+(4) + G+(3)G_(4)] (46a) 
J^l = G+(2)G_(3)G_(4) (46b) 
J'l ^ G+(2)G+(3)G+(4) (46c) 

Note that until the integration ranges for the momenta ^3? ^4 are specified, the same 
analytical expression F+(l) (^Sj) corresponds to the each of graphs depicted in Fig. 2. 

The following summation of Matsubara frequencies in (|45|) is straightforward. Due to 
conditions ( p6D we approximate the occurring hyperbolic functions up to exponentially small 
terms (9(exp(— A/T)). With this accuracy we drop all contributions from the term !F^. 
Analysis of the available kinematic configurations from Table 1 for other terms shows that 
some of them are also exponentially small, so we put 

(47a) 
(47b) 
(47c) 

Our calculations show that the term JPp* evaluated for the configurations of the graph a 
(solutions I, II) and c (I, II) gives four equal contributions to the self-energy. Analogously, 
two contributions of JF" from b (I, II) are equal to each other. In order to demonstrate how 
these calculations are done, let us consider as an example J^'^^ for the configuration (a, I). 
To obtain its total contribution to the self-energy, we have to to integrate two independent 
momenta within the ranges indicated in Table I for fixed n, m and then to sum over all 
available cascades (a, I) according to (P|a). This gives: ^ 

N-2m /,-A-mA+A pAn-mA+A ii 

^m) ^ = '^{£ + 3£)J:J: dkJ (48) 

«1 -"-O ^^-^ J-K-mA JA„-mA a+ — /v2 + /C4 



-res 
*^+ 


bI,II ~ 




aI,II ~ 




cI,II ~ 



"lower-order" cascades) are in fact recovered in the double integrals appearing in the limit A — > as some 
points or lines, i.e., manifolds of inferior measure. 
^ Strictly speaking, couplings in the sunrise graphs depend on the smallest shell number, i.e., n. By pulling 
them out of sums and integrals, as should be clear from the content of this section, we actually approximate 
them by their values at the largest shell number (i.e., at A). This approximation (called "local" in Ref.[^) 
is common for all RG schemes, and in any case possible corrections to it lie beyond the two-loop level. 
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where 



(49) 



In Eq. (ji^) and in the following we use dimensionless couplings absorbing the density of 
states factor in their definition, i.e., 



In the limit A — the sums above neatly combine into a simple double integral 



- !^ 

al 16 



which gives 



asi^(o) 



duj 



al 



-|(Ao+A) 

4- 16 



Ao+A+fc2 



dki 



k2 



(in^....) 



(50) 



(51) 



(52) 



where the omitted terms const + (9(A/Ao) are irrelevant for the RG flow equations. Com- 
bining evaluated in the same way other contributions and using the definitions (|40| ; PTD , we 
obtain the two-loop RG equation 



dl 



-^{gl + ?>gl + 2gl) 



(53) 



This equation was first obtained in Ref. 



C. Wilson-KadanofF vs field-theoretic schemes 

In order to get more insights from the obtained results, let us first rewrite some of the 
previous formulas. Using the transfer momentum vector g = ^2 — ^4 in Eq.(^), we get 

S?(l) = "-fM + 3g!) r dk r (54) 

In the above equation we dropped the index of the vector k2 and neglected terms A <^ Aq 
in the limits of integrals, since we are eventually interested only in the A — limit. Other 
three non-zero contributions [cf. Eq.(^7|a)] to the cs-term of the self-energy after analogous 



variable changes can be reduced to the form (|5^) but with different integration ranges. 



Combined together, these contributions can be written as: 

^ '^M + 3.3 *e(F > A^) r (55) 

where G is the event (Heaviside) function. 

Thus, by applying the canonical Wilson-Kadanoff scheme where all momenta in the 
diagrams lie strictly within an infinitesimal shell at each step of the RG procedure, and by 
summing all pertinent infinitesimal contributions through singling out kinematically-allowed 
cascades, we end up with results which could have been obtained in the field-theoretic 
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scheme. As one can see, formula (^) for the self-energy resembles the one derived in the 
RG scheme where momenta are let to run freely within a whole band restricted only by the 
IR and UV cutoffs. 

To sharpen the last statement, let us derive now the self-energy in the field-theoretic 
RG scheme with the cutoff regularization. (Thorough description of various field-theoretic 
RG approaches can be found, e.g., in Ref.[0].) Let us assume that we have the same bare 
effective action as we used before. Since the theory is logarithmical, we regularize divergent 
integrals by introducing IR and UV cutoffs. We assume that they satisfy the same conditions 
as those for the WK effective action (PBf). Since a priori the cutoffs in the field-theoretic and 



WK approaches are not the same, we will mark the former with tildes.^ Contrary to the WK 
scheme, in the field-theoretic approach the self-energy contribution we are looking for: (i) 
comes from the only one sunrise graph as those depicted in Fig. 2; (ii) independent momenta 
in the graph are restricted by the model cutoffs only ; (Hi) the graph has a combinatorial 
factor 1/2. All the summations before the final momentum integrations in this graph lead to 
the same formulas as in the previous subsection (except the overall factor 1/2). We choose a 
couple of independent momenta ki restricted by the cutoffs A, Aq. It can be shown that the 
term of the self-energy involving charge and spin couplings can be brought to the form (|55D 
with A I— > A, Aq 1-^ Aq up to 0{A/Aq). With the same accuracy the self-energy Umklapp 
terms coincide in the WK and field-theoretic approaches. Thus both schemes result in the 
same leading ln(Ao/A) term in the derivative of the self-energy with respect to the frequency, 
and therefore to the same RG equation (|53|). 

The relationship between the two RG schemes can be seen even in a simpler way in the 
limit A — i> 0, i.e., the fixed point limit we are always interested in. Let us draw the regions of 
momenta integrations corresponding to the solutions from Table I summed over available n 
and m. As a couple of independent momenta we choose (fcs, k^). These regions are shown in 
Fig. 3. Let us remind also that while showing kinematically non-equivalent cascades (a-c) in 
Fig. 2 we have used the (anti) symmetry of the vertex, and, correspondingly, that of graph's 
contributions under the exchange (3^4). So the sum of these three graphs can also be 
written as (a + b + c) = ^{a + b + c + a' + b' + c') , where the primed graphs correspond to the 
above exchange. As one can see from Fig. 3 the regions for unprimed and primed solutions 
(the primed solutions are obtained from their unprimed counterparts by the reflection over 
the main diagonal), fill completely the hexagon shown in bold. The latter is nothing but 
the whole phase space of the sunrise graph 

5restr : \ki\ <Ao, ^ = 2,3,4; fci = 0. (56) 

Now assuming that the summations and integrations commute we can write the total sunrise 
self-energy correction H^f^'^ resulting from the sum of three kinematically distinct graphs 2a-c 
over two solutions for each graph and over all available cascades as a single term 

j.casc(^) _ 1 /■ r(12,34)G(2)G'(3)G(4)r(34, 12) 5(1+1^(1 + 2-3-4) (57) 

2 J(2,3,4);fe,;g5rcstr 

Here we came back to the original thin variables comprising chiral and spin indices, Matsub- 
ara frequencies, and momenta {uji,ki 0). Thus the self-energy correction has the same 



^ For the IR cutoff A to be the smaUest scale of the problem we will work with the theory in the zero- 
temperature limit. For the following comparisons note that the second condition ( p6| ) effectively brings 
the WK RG results to the zero-temperature limit as well. 
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analytical form as the one written in the cutoff field-theoretic approach. In the latter how- 
ever, after removing one momentum by the delta-function (e.g., ^2) only two independent 
momenta are explicitly restricted by Aq. This results in the contributions from the extra 
states forbidden by the whole phase space constraints. These extra states lie in the triangles 
in the first and third quadrants of the {k^, k^j-plane (cf. Fig. 3). In the logarithmic theory 
as the one we are dealing with, integration over these forbidden states may result in only 
the finite terms to the self-energy correction, thus leaving the RG equation ( |5B| ) unaffected.^ 
The phase space mapping of Fig. 3 holds in the fixed-point limit A ^ at non-zero 
temperature as well. Thus the temperature effects can be calculated from the UV cutoff RG 
scheme upon integration of the all phase space restricted by the hexagon in Fig. 3.^ As it 
should be clear, this is not equivalent either to calculation within WK RG stopped at some 
/* when vfA{1*) ~ T, or within the two-cutoffs RG with the IR cutoff provided by T: such 
two ways of calculations give, strictly speaking, approximations only. 

IV. COUPLINGS AT THE TWO-LOOP LEVEL 

After explaining most of technicalities in the previous section, we will be brief in present- 
ing the results for the RG coupling equations. Let us first write the two-loop perturbative 
expansion for the four-point IPI vertex. It is depicted graphically in Fig. 4. (The tadpole 
self-energy corrections are absorbed in the one-particle Green's functions renormalizations.) 
We calculate the couplings from r(12,34) at 



where Uo is the minimal Matsubara fermionic frequency. The latter is set to zero after 
appropriate analytical continuation. The couplings are extracted from F by convolutions 
[cf. Eqs.(^) and below]. The limits (^) considerably reduce the number of cascades to 
take into account while calculating the renormalized couplings. 

One-loop corrections are generated by a simple contraction shown in the first line of 
Fig. lb. After summations over the shell number these contributions reduce to three one- 
loop graphs with the signs and factors as in the pertubrative expansion in Fig. 4 and the 
loop momentum A < \k\ < Aq. 

Pertinent double-loop graphs are generated by the 2-step cascade as follows: at the n-th 
step a simple one loop is generated (cf. first line of Fig. lb) and at the m-th step the left 
(or right) pair of legs of this graph is contracted with the four-leg vertex. The shell numbers 
n and m are independent and can take any value from 1 to N. The sums of these cascades 
result in the three double-loop graphs (DL) as they stand in the second line in Fig. 4 with 
the loop momenta restricted by A < \k\ < Aq. 



^ It is possible, (e.g., it happens for the component T,"^^) that summation over frequencies bringing the 
Fermi functions into play, suppresses completely the differences between two schemes coming from the 
extra states. For E'^'' the contributions from the first and third quadrants are completely eliminated 
by the Fermi functions [cf. Eq.(^7|a)]. But gencrically there are no reasons for that, and the Umklapp 
contribution provides a counterexample for E"^**. Cf. Eqs.(p7|b,c). 

^ Note, however that the explicit form of the flow equations, i.e., the Gell-Mann-Low functions, derived 
here in the limit T — > 0, are modified at T 7^ when decreasing cutoff A(Z) reaches the scale provided by 
the temperature. 




(58) 
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The two-loop graphs combined in the groups (ZZ, ZC, CZ) in Fig. 4 are generated by the 
3-step cascades. For concreteness let us take the first graph from the group ZZ. For each 
type of this loop configuration with a fixed position of external variables (top left in Fig. 5) 
there are six graphs corresponding to all possible types of the 3-step contractions shown 
in this figure. Direct inspection shows that the sign and combinatorial factor of all these 
six graphs are the same and coincide with those of their counterpart in the perturbative 
expansion in Fig. 4. Due to conditions (^) we retain only such cascades that the legs of the 
vertex a are contracted at the same step. From the conservation law the momenta k^^h^ 
are equal (opposite for the CZ graphs), so the phase space constraints near the vertex a 
are satisfied trivially. To find the allowed kinematic solutions from the second independent 
momentum conservation (e.g., at the vertex (5 of this graph) we can use directly the results 
of the previous section. Indeed, if we identify ^8)^5; ^7 of the graphs in Fig. 5 with 
ki,k2,k3, ^4 of Fig. 2, respectively, we recover the situation analysed previously for the self- 
energy graphs. The kinematic solutions for the graphs of Fig. 5 can be taken from Table I 
(in the limit L = 0) as follows: i i— > a, m t— > c'. Hi i— > a', iv i— > 6', f i— > c, vi b. Then 
contributions i — vi (cf. footnote on p. 0) can be combined into a single graph (top left in 
Fig. 5) with two independent momenta lying in the phase space shown in Fig. 3. Analogous 
conclusions can be reached for other graphs of the groups ZZ, ZC, CZ. Thus, provided the 
stipulated above phase space of each graph, the renormalization of the vertex after mode 
elimination between A and Aq (A 0) is given by Fig. 4. 

The rest of the calculations is straightforward. We introduce the renormalized (dimen- 
sionless) couplings as follows: 

g, = Z^Z^'gl'-^ (59) 
The vertex renormalization constants Z^-^ in terms of g^^ are found to be: 

Z-' = 1 + ^/ + kg! + 3g! - 2gl)l - ^3^(1 - (60a) 
9c 8 gi' 

Z:^ = l + gsl + \{gl - gl + 2gl)l (60b) 

Z-,' = \^gcl- \{9l - Sgl)l + \{gl - gl)l' (60c) 
Combining the equations 

i=c,s,3 ^ 



with Eqs.( |53| , |60D we obtain the RG flow equations for couplings: 



^9c 2 l 2 /CO N 

= 29c9-i (62a) 

^ = 9z9c - ^^3(^3 + 9I) (62c) 

These equations were obtained earlier within the field-theoretic RG analysis 0. Note that 
in the above equations we dropped all the contributions from g^ assuming that they can be 
absorbed in the renormalization of Fermi velocities. (See, e.g., Ref. and more references 
therein.) 
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V. CONCLUSIONS 



The equivalence of the RG flow equations in the WK and field-theoretic schemes at the 
two-loop level is proven explicitly for the ID g-ology model. To obtain the higher-order RG 
equations within the canonical WK approach, it is crucial to sum over all available cascades 
of contaractions occurring at the infinitesimal steps of the WK procedure. Moreover, in 
the fixed-point limit A — »• the diagrammatics of the WK RG summed over all cascades is 
shown to be the same as that of the field-theoretic RG with the cutoff regularization and 
restricted phase space. 

Generalizations for the cases of classical or quantum bosonic and fermionic fields and 
higher dimensions are straightforward within the developed technique, resulting in the same 
conclusion about the diagrammatics, provided that the relevant phase space of those fields 
(i.e., the low-energy sector) can be brought to the same origin. 

If the latter condition can not be satisfied easily, like, e.g, for the low-energy sector 
of higher dimensional condensed-matter fermions possessing a Fermi surface, the proof of 
equivalence is not a straightforward extension within the proposed approach, and more 
explicit work is needed to demonstrate it. However we conjecture that even in that case 
the summation over cascades within the WK RG tantamount to applying the field-theoretic 
RG with a (bare) UV cutoff for the low-energy sector and locally-imposed phase-space 
constraints of the integrated momenta. (For the examples of how these constraints may be 
realized, see the review by Shankar 0.) 

At the level of three loops and higher, different RG schemes may result in different flow 
equations. Within the present approach, it follows from the necessity to take into account 
the dependence of running couplings on the shell number in higher orders (cf. footnote 
on p. |1^). So various cascade contributions could not be "compactified" in a single "field- 
theoretic" graph as it occurs at the two-loop level. The latter conclusion is hardly surprising, 
since the differences in the Gell-Mann-Low functions in higher orders are known to occur 
even within different versions of the field-theoretic RG. (See, e.g, Ref.|]14|. In the context 



of the ID (yf-ology, see also Ref.[T^ for the discussion of the bandwidth- and transfer-cutoff 
regularizations in higher orders of RG.) 

In condensed-matter problems, where the bare UV cutoff Aq might be a meaningful 
scale and other scales (e.g., T, Kp) may be present as well, the WK RG scheme should be 
applied, as the one which is just an exact way to to integrate the partition function [0, 
The WK RG does not rely on any assumptions of the scales involved, scaling/logarithmicity 
of theory, etc, and includes contributions only from the states lying within the phase space 
of the problem. This allows to control the results (at least in principle) in all orders. In 
practical terms, after summations over all cascades the WK scheme asymptotically maps 
to the much easier to handle field-theoretic RG (up to three loops) with forbidden states 
excluded. We did not verify whether possible simplifications of the WK scheme can be done 
in even higher orders, but going beyond two-loop level is a formidable problem in any case. 
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(a) 




n+m 



(b) 

FIG. 1: Diagrammatic representation of the two types of cascades of contractions in the second 
order over interaction, (a): 3-step cascade; contractions are done at the steps n,n + m,n + m + k. 
(b): 2-step cascade; contractions are done at the steps n and n + m. 




(d) (e) 



FIG. 2: Two-loop kinematically-distinct self-energy diagrams generated after TV steps of the 
Wilson-Kadanoff procedure. Graphs (a — c) are generated by the 3-step cascade (see Fig. la), 
n < n + m < N — m'; (d-e) by the 2-step cascade (see Fig. lb). Integrcs shown near the composite 
variables (1 — 4) indicate the shell number where the corresponding momenta are restricted to lie. 
The signs and combinatorial factors of the graphs stand as they enter the action's binomial term 
?/'(l)'0(l). Implicit sum over all allowed n,m{m') is assumed. 
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-Ac 

FIG. 3: Regions of momenta integrations corresponding to the solutions from Table I in the hmit 
A ^ 0. The regions for the primed solutions (not given explicitly in Table I) are obtained by 
reflection over the main diagonal (bold dashed line), so the whole set of 12 solutions fill completely 
the restricted phase space (|56D, i.e., the hexagon shown by the solid bold line. The vertical thin 
dashed lines are guides for an eye. 



(ZZ) 



(ZC) 



(CZ) 




FIG. 4: Two- loop perturbative expansion for the four-point vertex. The tadpole self-energy 
corrections are absorbed by the Green's function renormalization. The vertex variables indicated by 
numbers comprise the complete set: chiral and spin indices, Matsubara frequencies, and momenta. 
For convenience the two-loop graphs are separated in four groups (DL, ZZ, ZC, CZ) discussed in 
the text. Note that in the graph CZ we depicted a twisted particle-hole loop, not a double loop. 
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FIG. 5: 3-step cascades of contactions generating two-loop corrections (topology of the ZZ, ZC, 
CZ graphs) to r(12, 34) in the limit ([58|). As an example, six distinct cascades corresponding to the 
kinematics and the external legs configuration of the first graph from the group ZZ of Fig. 4, are 
shown. For the sake of space, the vertex labels and notations for the integrations dummy variables 
discussed in the text (the same for all depicted graphs), are explicitly indicated in the top left one 
only. The integers in the graphs (i-vi) show the shell numbers where momenta (fcs, .., /cg) lie. Note 
that as in the case of the self-energy corrections (cf. Figs. 1-2), the 3-step cascades come from 
higher many-body interactions, themselves generated along the RG flow. In the graph (i), taken as 
an example, the first contraction (step n) leads to a four-particle interaction, the second step n+m 
to a three-particle one, and the final third step N-m yields the two-particle interaction correction 
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